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ABSTRACT 



Context. The recent precise measurement of the mass of pulsar PSR Jl 6 14-2230, as well as observational indications of even more 
massive neutron stars, has revived the question of the composition of matter at the high densities prevailing inside neutron-star cores. 
Aims. We study the impact on the maximum possible neutron-star mass of an "exotic" core consisting of non-nucleonic matter. For 
this purpose, we study the occurrence of a first-order phase transition in nucleonic matter. 

Methods. Given the current lack of knowledge of non-nucleonic matter, we consider the stiffest possible equation of state subject only 
to the constraints of causality and thermodynamic stability. The case of a hadron-quark phase transition is discussed separately. The 
purely nucleonic matter is described using a set of unified equations of state that have been recently developed to permit a consistent 
treatment of both homogeneous and inhomogeneous phases. We then compute the mass-radius relation of cold nonaccreting neutron 
stars with and without exotic cores from the Tolman-Oppenheimer-Volkoff equations. 

Results. We find that even if there is a significant softening of the equation of state associated with the actual transition to an exotic 
phase, there can still be a stiffening at higher densities closer to the center of the star that is sufficient to increase the maximum possible 
mass. However with quarks the maximum neutron-star mass is always reduced by assuming that the sound speed is limited by c/ V3 
as suggested by QCD calculations. In particular, by invoking such a phase transition, it becomes possible to support PSR J1614- 2230 
with a nucleonic equation of state that is soft enough to be compatible with the kaon and pion production in heavy-ion collisions. 

Key words. Stars: neutron - equation of state - dense matter - gravitation - methods: numerical 
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1. Introduction 

Born from the catastrophic gravitational core collapse of mas- 
sive stars (mass M > 8M ) at the end point of their evolution, 
neutro n stars are among th e most compact objects in the uni- 
verse (lHaensel et all I2007I) . with a radius of the order of 10 
km and a mass of around 2 M . A few meters below the sur- 
face at densities above ~ 10 4 g.crrT 3 , matter is so compressed 
that all atoms are fully ionized, and are arranged on a regular 
Coulomb lattice of nuclei, neutralized by a gas of degenerate 
electrons; this is the "outer crust". Deeper in the star, nuclei be- 
come more and more neutron rich, as a result of electron cap- 
ture, and at a density of about ~ 4 x 10 11 g.cnT 3 neutrons be- 
gin to drip out of the nuclei. This marks the transition to the 
inner crust, an inhomogeneous assembly of neutron-proton clus- 
ters and unbound neutrons, neutralized by the degenerate elec - 
tron gas (Pethi ck & RavenhaUlll995tlChamel & HaenselL [20081. 
When the density reaches about ~ 10 14 g.crrT 3 (about half the 
density found at the center of heavy nuclei), the crust dissolves 
into a uniform plasma of neutrons with a small admixture of pro- 
tons, neutralized by electrons and, at slightly higher densities, 
muons. 

The mass and radius of a neutron star are determined by the 
equation of state (EoS) over the full range of densities found 
in the star, i.e., by the relation between the pressure P and the 
mass-energy density p, although the core will play a dominat- 
ing role. This has motivated many studies of purely nucleonic 
neutron-star matter (N*M), i.e., a homogeneous and electrically 
charged neutral liquid of nucleons and leptons in beta equilib- 



rium. These studies consist of simple extensions of the large 
number of many-body calculations performed since the begin- 
ning of the 1950s on so-called nuclear matter, consisting of 
just neutrons and protons that interact via "realistic" forces fit- 
ted directly to experimental nucleon-nucleon phase shifts and 
to the properties of bound two- and three-nucleon systems (the 
Coulomb force being switched off). The EoS of purely nucleonic 
N*M has been determined in such many-body calculations up to 
the highest densities found in the most massive neutron stars. 
The maximum neutron-star mass obtained using different many- 
body methods and realis tic forces is predicted to lie in the range 
between -1.8 -2.5M Q (iLi & Schu lze. 2008; FuchsL l2008h and is 
therefore compatible with the measured value 1 .91 + 0.04M o fo r 
the mass of the pulsar PSR J1614-2230 dDemorest et all l2010). 

However, the core of massive neutron stars is likely 
to contain not only nucleons and leptons but also other 
particles like hy perons, meson condensates, or even de - 
confined quarks ilPage & Reddvl l2006t IWeber et all 120071) . 
Neutron stars with a hyperon core are s ometimes referred 
to as "hyperon stars" dGlendenningl l2000l) . Now according 
to Brueckner-Hartree- Fock calcula t ions using realistic two- 
and three-body forces (|Vidana et al.l 1201 U iBurgio et al.l 1201 U 
Schulze & RiikenL 1201 ll) . the appearance of hyperons in dense 
matter softens the EoS considerably thus lowering the maximum 
neutron-star mass to an almost unique value around 1.3 — 1 .4M . 
On the other hand, some relativistic mean-field calculations in- 
cluding hyperons ca n support neutron stars that are as massive as 
PSR J1614-2230 dBednarek et aUl2012tlSulaksono & Agrawa " 
1201 2t iJiang et all l2012t IWeissenborn et all 120121; IZhao & 
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120121) . This discrepancy can be understood at least partly from 
the fact that the maximum mass is very sensitive to the vari- 
ous hyperonic couplings, and these are determined very poorly 
since the limited nuclear and hypernuclear data that are relevant 
constrain the EoS only in the vicinity of the saturation density, 
whereas the maximum neutron-star mass is mostly determined 
by the EoS at much higher densities. Likewise, some relativistic 
mean-field models including meson conde nsates are able to pre- 
dict th e existence of massive neutron stars (Gupta & Arumugam, 
120121) . 

Another possibility for raising the maximum mass of neu- 
tron stars above ~ 2M lies in the deconfinement of quarks 
in t he core (such stars a re generally called "hybrid stars", see 
e.g. iGlendenningl [2000). However, most of the calculations on 
which this conclusion is based are phenomenological in the 
sense that they lack a direct r elationship w i th qua ntum chro- 
modynamics (QCD): see, e.g. Alfor d etldl d2007l) and refer- 
ences therein. Even though recent perturbative QCD calcula- 
tions lead to predicting compact stars compatible with PSR 
J1614- 2230 dKurkela et all 1201 Oh . these calculations are not 
strictly valid for the densities prevailing in neutron stars, even 
in the most massive ones. In any case, whether the densities 
reached in neutron stars are high enough for deconfinement to 
occur is still an open question. 

The EoS of neutron-star cores allowing for the presence of 
all kinds of particles (nucleons, leptons, hyperons, meson con- 
densates, deconfined quarks, etc.) thus remains highly uncer- 
tain, and this situation is unlikely to be changed in the near fu- 
ture. Indeed, understanding the properties of high-density mat- 
ter would require a consistent treatment of the various hadron 
species taking their internal structure into account based on 
QCD. Unfortunately, solving the equations of QCD in the non- 
perturbative regime prevailing in neutron-star cores appears as 
an extremely challenging problem. In comparison, the ambi- 
guities in the EoS of purely nucleonic matter are much less 
acute, despite the divergences between different calculations at 
the higher densities found in neutron-star cores. 

In this paper we do not deal with all the large uncertainties in 
the underlying physics of the EoS of the "exotic" non-nucleonic 
matter that might be found in neutron-star cores. Rather, we 
determine in all generality the optimal possible increase in the 
maximum neutron-star mass over what is found with purely nu- 
cleonic N*M. We only impose the constraints of causality and 
thermodynamic stability, i.e. the condition that at a given pres- 
sure a phase transition will occur only if the Gibbs free energy 
per nucleon is lowered. For the purely nucleonic N*M, we use a 
set of unified EoSs based on the n uclear energy-density func- 
tional theory (Go rielv et all 1201 Ol) . as described in Sectj2] In 
Sect[3] we discuss the thermodynamics of a possible transition 
to an exotic phase, and then the EoS of such a phase. For the 
latter we consider first the case where the stiffness of the ex- 
otic phase is only limited by causality, i.e., the requirement that 
the speed of sound cannot exceed the speed of light, c. We then 
take account of the fact that, according t o both perturba tive QCD 
calculations at zero temperature (Kur kela et all [2010) and non- 
pertu bative lattice QCD calculations at finite temperatures (see 
e.g. lKarscRl2007l : lBorsanvi et all l2010. and references therein), 
the speed of sound in a gas of deconfined quarks cannot exceed 
c/ V3, and accordingly modify our causally limited EoS. The 
implications for the maximum neutron-star mass are discussed 
in Sect|4] while in Sect|5]we summarize our conclusions. 



2. The nucleonic equation of state 

To assess the role of a transition to a non-nucleonic phase we 
must begin with a purely nucleonic EoS that has been well 
adapted to the description of neutron stars whose cores are as- 
sumed to be non-exotic. Suitable such starting points are pro- 
vided by the family of three EoSs that we have developed 
to provide a unified treatment of all parts of neutron stars. 
These EoSs are based on nuclear energy-density functionals 
that have all been derived from effective interactions that are 
generalizations of the conventional Skyrme forces in that they 
contain terms that depend s imultaneously on momentum and 
density (iGorielv et all 12010). The parameters of this form of 
Skyrme force were determined primarily by fitting measured 
nuclear masses, which were calculated with the Hartree-Fock- 
Bogoliubov (HFB) method. For this it was necessary to sup- 
plement the Skyrme forces with a microscopic contact pairing 
force, phenomenological Wigner terms, and correction terms 
for the spurious collective energy. In fitting the mass data, we 
simultaneously constrained the Skyrme force to fit the zero- 
temperature EoS of homogeneous pure neutron matter (NeuM), 
as determined by many-body calculations with realistic two- and 
three-nucleon forces. Actually, several such calculations of the 
EoS of NeuM have been made, and while they all agree very 
closely at nuclear and subnuclear densities, they differ in the 
stiffness that they predict at the much higher densities that can 
be encountered towards the center of neutron stars. Functional 
BSkl9 was fitted to a soft EoS of NeuM (the one labeled "UV14 
plus TNI" in IWiringaetal.il 19881 combined with the EoS of 
Friedman & Pandharipande 1981) , BSk21 to a very stiff EoS 
(the one labeled "VI 8" in lLi & Schulzd2008l) . while BSk20 was 
fitted to an Eo S of intermediate st iffness (the one labeled "A18 + 
5 v + UIX*" in lAkmal etai1ll998l which we abbreviate as APR). 
All thre e EoSs are con s istent with quantum Monte Carlo calcu- 
lations (Gando lfi et all 120121) . Even though our EoSs were not 
fitted to symme tric nuclear matter they a re all consistent with 
the constraint of Daniel ewicz et all d2002l) deduced from heavy- 
ion collisions. Furthermore, the strength of the pairing force at 
each point in the nucleus in question was determined so as to ex- 
actly reproduce realistic 1 Sq pairing gaps of homogeneous nu- 
clear ma t ter of the appropriate density and charge asymmetry 
(IChamelL l2010l) . Finally, we imposed on these forces a number 
of supplementary realistic constraints, the most notable of which 
is the suppression of an unphysical transition to a spin-polarized 
configuration, both at zero and finite temp eratures, at densities 
found in neutron stars and supernova cores dchamel et ail [2009; 
IChamel & GorielvLl2010HGorielv et all l2010h . The form of our 
functionals was flexible enough for us to satisfy all these con- 
straints and at the same time fit the 2149 measured masses of nu- 
clei with N and Z > 8 g i ven in the 2003 Atomic Mass Evaluation 
(AME) (WaDs traetal.U 2003) with an rms deviation as low as 
0.58 MeV for all three models, i.e., for all three options for the 
high-density behavior of NeuM. 

These functionals are very well adapted to a unified treat- 
ment of all parts of purely nucleonic neutron stars, given not 
only the NeuM constraints to which they have been subjected 
but also the precision fit to masses, which means that the pres- 
ence of inhomogeneities and of protons is well repr esented. We 
used these functionals in lPearson et al .1 (120111120121) to calculate 
th e properties of the ou ter and inner crusts, respectively, while 
in Cha mel et all (1201 lh we determined the maximum possible 
neutron-star mass for each functional. For this calculation we 
had to solve the Tolman-Oppenheimer-Volkoff (TOV) equations 
(lTolmanlll939HOppenheimer & Volkofail939l) for different val- 
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Fig. 1. Left panel (a): range of energies per baryon (defined by 
e = S/n - m n c 2 where m„ is the neutron mass) of nucleonic 
matter in beta equilibrium at zero temperature as a function 
of the baryon densi ty for the unified Brussels-Montreal EoSs 
(lGorielvetal.Ll20Tol) . The shaded area reflects the different de- 
grees of stiffness of these EoSs: the lower limit corresponds 
to BSkl9, the upper to BSk21. Right panel (b): corresponding 
range of pressures. 
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Fig. 2. Range of neutron -star masses and rad ii for the unified 
Brussels-Montreal EoSs (lGorielvetal.Ll2010h . The shaded area 
reflects the different degrees of stiffness of these EoSs: the lower 
limit corresponds to BSkl9, the upper to BSk21. For compari- 
son, we have indicated the meas ured mass of PSR J16 14-2230 
including estimated errors from iDemorest et al.l (t2010b . 



ues of the central density, thereby obtaining the mass M of the 
star as a function of its radius R. The total mass of a neutron star 
depends largely on the core properti es (although we always ac- 
count for the crust using the EoS of iPearson et al.l I20TTII20T2L 
for the appropriate functionals) and thus on the EoS of homo- 
geneous N*M: Fig. Q] shows the range of uncertainty spanned 
by our three functionals, BSkl9 being the softest and BSk21 the 
stiffest. The corresponding uncertainties in the neutron-star mass 
for a given radius are shown in Fig. [2] from this same figure we 
can infer that BSk21 is stiff enough at high densities to support 
neutron stars as massive as PSR J1614-2230, but that functional 
BSkl9 is too soft (BSk20 can also support PSR J1614-2230). It 
is worth notin g that our functional s are consistent with the radius 



3. Transition to non-nucleonic phase in dense 
matter 

3. 1 . Thermodynamic equilibrium 

In view of the current uncertainties about the composition of 
matter at supernuclear densities, we simply assume that above 
a baryon density n^, nucleonic matter undergoes a first-order 
phase transition to some unknown exotic phase with a baryon 
density «x > «n- The exotic phase could consist of a pion 
condensate, a ka on condensate , hyperonic matter, or deconfi ned 
quarks (see e.g. iGlendennTn g, 2000; Hae nsel et al.L I2007I) . In 
the region of coexisting phases, «n < n < rix, thermodynamic 
equilibrium requires the constancy of the pressure P and baryon 
chemical potential jj. : 



Stein er et al l (§010) inferred from observations of p <=xo(«) - JVcOn), Mexo(«) - /w("n) . 



constraints of 
X-ray bursters and low-mass X-ray binaries. 

Since the softest of our Skyrme functionals that is com- 
patible with the mea sured mass of PSR J1614-2230 is BSk20 
dChamel et all 1201 lh . the question arises as to whether it is 
not too stiff to be compa t ible with the analysis of K + pro- 
ductio n dFuchs et all l2Q0lb ISturm et al~L 1200 t iHartnack et ail 
2006) and 7r~/7r + production ratio (IXiao et al.L 120091) that have 
been measured in heavy-ion collisions. In particular, the former 
analysis suggests that the EoS of symmetric nuclear matter is 
much softer than what is obtained with BSk20 over the range 
2«o ^ n < 3no, whereas the latter analysis concludes that over 
the range 2«o 55 n < 3.5«o the symmetry energy rises signif- 
icantly less steeply than predicted for the APR EoS, on which 
BSk20 is based. These results, taken at face value, discriminate 
against both BSk20 and BSk21, and favo r BSkl9. Various ex- 
otic mechanisms, such as a "fifth force' ' (IWen et all [2 009) or 
variations of the gravitational constant (IWen et all 12012b . have 
been proposed to simultaneously account for both this result 
and the existence of PSR J1614-2230, but we shall see that the 
much more prosaic explanation of a transition to a non-nucleonic 
phase may suffice. 



(1) 



where the subscripts "exo" and "nuc" are used to denote the ex- 
otic and nuclear-matter EoS, respectively. As a result, using the 
general expression 



8 = n/j.-P 



(2) 



valid at the zero temperature that we assume throughout this pa- 
per, we see that the energy density of the coexisting phases varies 
linearly with the baryon density in the range < n < nx ac- 
cording to 



£exo(n) = «j"nuc(«N) ~ ^nuc(«N) • 



(3) 



At the density «x, which we always suppose to be lower than the 
central density « cen of the star, the phase transition has reduced 
the pressure from -P n uc(«x) to P nuc ("N)- 

In this simple picture of a first-order phase transition, the 
two phases cannot coexist in the star because the densest phase 
will sink below the other. In hydrostatic equilibrium, the two 
phases would thus be spatially separated with the density vary- 
ing discontinuously at the interface. In reality, the two distinct 
phases will generally rearrange themselves by forming a mixed 
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phase in a finite region of the star (Tatsumi et al.L I20TTT) unless 
the surface tension and the Coulomb interaction are sufficiently 
strong, a s could be the case for the hadron-quark phase transition 
(see e .g. lHeiselberg et al.Lll993tlAlford et al.Ll200ll:lEndo et all 
2006). The presence of a mixed phase leads to a smooth tran- 
sition with the pressure increasing monotonically with the den- 
sity instead of remaining constant. This situatio n generally in- 
creas es the maximum neutron-star mass (see e.g. iGlendenningl 
2000, section 9.3). However, given the uncertainties pertaining 
to the EoS of the mixed phase, here we suppose the least fa- 
vorable case of a first-order transition considered in Eqs. ([TJ 
and (f3]l w hich will provide a lower bound on the maximum mass 
(IRhoades & RuffiniLll974l) . 

3.2. Equation of state for exotic matter 

Causal limit of EoS. Instead of considering specific models 
of non-nucleonic matter, we suppose that the EoS of the ex- 
otic phase is the stiffest possible, being only constrained by Le 
Chatelier's principle 



dP 



> 



and the causality requirement 
— < 



(4) 



(5) 



Since the maximally com pact stars are obtained for dP/dS - 1 
(IRhoades & RuffiniL 11974 . we will assume that the pressure of 
the exotic phase at densities n > rix can be expressed as 

.Pexo(«) = <5exo(«) ~ <5 e xo(«x) + J»uc(»n) . (6) 

recalling that P„uc(«n) = Pexo(«x)- After integrating the pressure 
by using Eq. (f2]i with fi = dS/dn, we find for the energy density 
of the exotic phase 



<5exo(«) 



«x|W(«n) 



1 + 



"X 



•Pnuc(«N) ■ 



(7) 



Using Eqs. (fj), ©, and (fT]), the baryon chemical potential in the 
exotic phase can be expressed as 



/"exoCP) =/W("n) a 1 + 



2(P-P nuc (« N )) 



«x/w("n) 



(8) 



From both Eqs. (0 and © it is seen that the EoS of the exotic 
phase is completely determined by the parameters «n and rix, 
and for any given pair of values of these parameters we have to 
determine the maximum possible neutron-star mass. But in this 
very general study we have a high degree of freedom in choosing 
the values of these parameters, and we do so in such a way as to 
obtain a maximum maximorum in the neutron-star mass, subject 
to certain physical constraints that we now discuss. 

Dealing first with n^, we find that increasing values of this 
parameter lead to decreasing values of the maximum mass, as is 
to be expected intuitively. We therefore consider the lowest pos- 
sible density «n consistent with nuclear data. In particular, we 
set «n = 0.2 fin -3 , which is slightly higher than the highest den- 
sity found in nuclei, as predicted by HFB calculations on more 
than 8000 nuclefl. Since much higher densities can be reached 



in heavy-ion collisions, the lack of any evidence of phase transi- 
tions might suggest that «n » 0.2 firT 3 . In fact, this is not neces- 
sarily the case because the conditions prevailing in neutron-star 
interiors are very different from those encountered in heavy-ion 
collisions. We do not discuss here the possibility of strange stars 
for which the exotic phase (strange-quark matter in this case) 
would be p resent in the entire star. At d ensities n < n^, we use 
the EoS of iPearson et al.l d201 ll 1201 2|) for the outer an d inner 
parts of the crust and the EoS of iGorielv et al.l (1201 Ol) for the 
purely nucleonic part of the core. 

As for «x, we note first that the pressure of the exotic phase 
remains lower than that of the nucleonic phase from the baryon 
density «n up to some density rip > rix, but that it is higher there- 
after, as shown in Fig. [3] Therefore the impact of a phase transi- 
tion will be to increase the maximum neutron-star mass provided 
the central density n cen is substantially higher than rip. However, 
rip cannot be freely adjusted because of the requirement that the 
exotic phase should be energetically favored. More precisely, at 
a given pressure P, the equilibrium phase is found by minimiz- 
ing the Gibbs free energy per baryon, which coincides with the 
baryon chemical potential. Now the baryon chemical potential 
of the exotic phase will rise more steeply than that of the nucle- 
onic phase, and at some pressure Pc the two phases will have the 
same baryon chemical potentials, // e xo(Pc) = /w(Pc)- This will 
occur when the density of the exotic phase reaches the value nc, 
such that P eX o(«c) = Pc- F° r pressures P > Pc, or equivalently 
for densities n > nc, the ground state of matter will once again 
be purely nucleonic. In order to exclude this unlikely possibility, 
we must have Pc higher than the central pressure P cen (or equiv- 
alently «c higher than n cen ) in the most massive neutron stars. 
Fixing Pc then completely determines the density rix, which is 
given by 



«x 



2(P C - Pnuc(» N )) 



/w(«n) 



A<nuc(Pc) 

/w(«n) 



1 



(9) 



where we have used Eq. ([§}. At the same time, we confirmed 
the intuitively plausible result that the lower Pc is, the greater 
the maximum possible mass, and accordingly we arranged for 
Pc to be very close to P cen . Now the maxi mum central pressure 
in a n eutron star is approximately given by (Lattimer & Prakash, 
l20TD) 

P cen « 2.034fi exo (« x ) = 2.034 (hx/w(»n) - A,uc(«n)) , (10) 

where we have used Eq. ((TJ. The optimum value of rix was there- 
fore obtained by solving Eqs. <j9j and (ITOl l with P cen = Pc- 
However, since Eq. ( fTQb is only approximately valid, we sub- 
sequently refined the value of Pc by solving the TOV equa- 
tions numerically in order to obtain a more accurate estimate 
of P C en- With Pc and nx determined, the density nc (for which 
Pexo(«c) = Pc) can be obtained from Eqs. © and (0. The cor- 
responding EoSs are as shown in Figs. [3] and |4] In particular, 
we note that at a given pressure the phase transition leads to a 
lowering of the baryon chemical potential. 

The density rip is completely determined by the nucleonic 
EoS, «n and ric, as follows. We first note that for the EoS given 
by Eq. © we have // = dS/dn - dP/dn. Differentiating Eq. (f2]i 
thus leads to 
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d[i dn 
/j n 

which after integration using Eq. ([T), yields 

(«p) jUexo(nx) («n) 



«x 



Tlx 



(ID 



(12) 
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Fig. 3. Upper panels: pressure as a function of the baryon num- 
ber density (up to the highest central density n cen found in the 
most massive neutron stars) for the softest (BSkl9, panel a) and 
the stiffest (BSk2 1 , panel b) of our nucleonic EoSs (dashed lines) 
and for the corresponding EoSs with a causally limited phase 
transition (solid lines). Lower panels: energy per baryon (defined 
by e — S/n - tn n c 2 where tn n is the neutron mass) as a function 
of the baryon density for the softest (BSkl9, panel c) and the 
stiffest (BSk21, panel d) of our nucleonic EoSs (dashed lines) 
and for the corresponding EoSs with a causally limited phase 
transition (solid lines). The shaded areas indicate the region of 
phase coexistence at densities «n < « < «x- 
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Fig. 4. Pressure as a function of the baryon chemical potential 
for the softest (BSkl9, panel a) and the stiffest (BSk21, panel b) 
of our nucleonic EoSs (dashed lines) and for the corresponding 
EoSs with a causally-limited phase transition (solid lines). The 
horizontal dotted line indicates the central pressure in the most 
massive neutron stars. 



Case of quark matter. The foregoing analysis may be invalid 
and the limiting EoS much softer if quark deconfinement takes 
place in the core of a neutron star. Accordi ng to both perturba- 
tive QCD calculations at zero temperature (Kurk ela et al.L [2010) 
and n onperturbative lattice QCD calculatio ns at finite tempera- 
tures dKarschl 120071 : iBor sanvi et all 1201 Ol) . the speed of sound 
in quark matter is limited by c/ V3. Assuming that this result re- 



mains valid in the interior of neutron stars, the stiffest possible 
EoS is given by 

•PquarkO) = T [SquarkO) ~ SquarkOx)] + Aiuc(«n) ■ (13) 



After integrating the pressure, we find for the energy density 



£quark(") = ^«X/Unuc(«N) 



4/3 



Aiuc(«n) , 



(14) 



instead of Eq. (jTj. This equation of state ressembles that ob- 
tained with the simplest MIT bag model, which has been com- 
monly used for describing the interior of compact stars (see e.g. 
iGlendenningl l2000t lHaensel et all l2007h . In this model, quarks 
are treated as massless and noninteracting particles confined in- 
side a "bag". The pressure of the quarks is then given by 



.Pquark(«) = -(SquaikO) ~ 45 j , 



(15) 



where B is the bag pressure. Comparing Eq.rtTJt and (TT5t yields 
the effective bag pressure 



B = -«x/W(«n) - Aiuc(«n) ■ 



(16) 



As a matter of fact, Eq.(14) is also a fair ly good approxima- 
tion of more realistic quark m atter EoSs (Zdu nik et all [2000; 
iGondek-Rosiriska et all l2000h . Using Eqs. dT3} and (El), the 
chemical potential is given by 



^quark(-P) = /W("n) 



MP - A,uc(«n)) 



jUnuc(«N)«X 



+ 1 



1/4 



(17) 



Again, the EoS of the quark phase is completely determined 
by the parameters «n and tlx, whose values must be specified. 
We choose the former as before, remarkin g that according to 
percolation simulations (Magas & SatzL 120031) the density «n at 
which the hadron-quark phase transition begins could be very 
close to the saturation density. The density nx is still determined 
by maximizing the neutron-star mass, subject to the constraint 
that there be no reconversion of the exotic phase into nucleonic 
matter at the highest densities prevailing in the most massive 
neutron stars. However, the situation is now more complicated 
than in the previous case, where the stiffness of the EoS is limited 
only by causality. The EoS for quark matter is shown in Figs. |5] 
and|6] where it will be seen that it is possible to have equality of 
both the chemical potentials and the pressures of the two phases 
at the same density tic ■ 



P quark(«c) = -P n uc(«c) . Mquark("c) = ^nuc("c) • 



(18) 



It follows from Eq. (j2]i that the energy densities of the two phases 
are also equal 



£quark(«c) = £nuc(«c) • 



(19) 



We found numerically that this situation corresponds to the max- 
imum possible neutron-star mass. It particularly needs to be em- 
phasized that while tic is now less than « cen , the exotic phase is 
not reconverted into nucleonic matter: although it is on the point 
of doing so at tic, the EoS of quark matter immediately soft- 
ens as n increases beyond tic- In any case, we see that while the 
pressure in the quark phase is higher than in the nucleonic phase 
at relatively low densities, the reverse is the case at the higher 
densities relevant to the maximum neutron-star mass that can be 
supported. We may thus anticipate that quark deconfinement will 
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Fig. 5. Upper panels: pressure as a function of the baryon num- 
ber density (up to the highest central density n cen found in the 
most massive neutron stars) for the softest (BSkl9, panel a) 
and the stiffest (BSk21, panel b) of our nucleonic EoSs (dashed 
lines) and for the corresponding EoSs with a quark phase tran- 
sition (solid lines). Lower panels: energy per baryon (defined by 
e = S/n - m n c 2 where m„ is the neutron mass) as a function 
of the baryon density for the softest (BSkl9, panel c) and the 
stiffest (BSk21, panel d) of our nucleonic EoSs (dashed lines) 
and for the corresponding EoSs with a quark phase transition 
(solid lines). The shaded areas indicate the region of phase co- 
existence at densities «n < « < «x- 



1000 



> 
a, 




1000 



1500 2000 

H [MeV] 



1000 



1500 2000 

H [MeV] 



Fig. 6. Pressure as a function of the baryon chemical potential 
for the softest (BSkl9, panel a) and the stiffest (BSk21, panel b) 
of our nucleonic EoSs (dashed lines) and for the corresponding 
EoSs with a quark phase transition (solid lines). The horizon- 
tal dotted line indicates the central pressure in the most massive 
neutron stars. 

Table 1. Stability of nonrotating neutron stars with exotic cores: 
range of central pressures of stars that are unstable with respect 
to radial oscillations, corresponding range of masses and radii. 



reduce the maximum neutron-star mass, assuming that the speed 
of sound is limited by cj y3. 

Concerning the low-density phase transition from nucleonic 
to quark matter, we see on comparing Figs. [3]and[5]that the den- 
sity range «x - "n over which the two phases coexist is much 
narrower than in the case of the transition to the causally limited 
EoS. 

The effective bag constants associated with the nucleonic 
EoSs BSkl9, BSk20 and BSk21 (respectively 78.6, 65.5, and 
56.7 MeV frrT 3 ) lie in the range of values that have been gen- 
erally adopted for studies of hybrid stars (see e.g. lHaensel et all 
2007). 



4. Neutron-star maximum mass 

For each of our three unified EoSs we have solved the TOV equa- 
tions that des cribe the global structure of spherical nonrot ating 
neutron stars (Tolman, 1939; Qppenhei mer & Volkoffl[T939h . As 
shown in a previous paper (IChamel et all 1201 ll) . the impact of 
the rotation on the maximum neutron star mass is negligibly 
small for stars having rotation periods comparable to that of PSR 
J1614-2230. For simplicity, we have therefore ignored rotation. 
The resulting mass-radius relations are plotted in Figs. |7]and[8] 
Neutron stars with central densities « cen > «x may be unsta- 
ble with respect to radial oscillations due to the strong soften- 
ing accompanying the phase transition in dense matter. Unstable 
confi gurations, which ar e characterized by the inequality (see 
e.g. lHaensel etalll2007l) 



'cen 



Force 


P cen (MeV fnr J ) M/M Q R (km) 


BSkl9 
BSk20 
BSk21 


3.32- 11.0 0.145 -0.188 14.5 - 15.5 
4.17- 9.28 0.207 -0.233 13.1 - 13.9 
4.87- 9.49 0.265 -0.290 12.6- 13.2 



Table 2. Same as Table[T]for neutron stars with quark cores. 



Force 


P cc „ (MeV fnr J ) M/M Q R (km) 


BSkl9 


3.32- 11.0 0.186-0.188 15.32- 15.33 



are indicated in Tables Q] and |2]for each of our three functionals. 
For neutron stars with quark cores, no instabilities were found 
for BSk20 and BSk21, and unstable configurations for BSkl9 
were found to be restricted to a very narrow range of masses and 
radii. 

The numerical values of the maximum neutron-star masses 
M max are indicated in Tables[3]and[4]for each of our three func- 
tionals. We also show in these tables the corresponding radius, 
the highest baryonic density «n of the nucleonic phase, the low- 
est baryonic density «x of the exotic phase, and the central bary- 
onic density « cen . 

In Table [5] we summarize the results of Cha mel et al 1 (120111) 
for the case of no exotic phase. Comparing with Table [3] we 
see that the result of allowing a transition to an exotic phase 
whose EoS is limited only by causality is always to increase 
the maximum possible neutron-star mass. Of particular interest 
is the case of functional BSkl9, which can now support pulsar 
PSR J1614-2230. Thus this functional allows us t o reconcile 
the existence of thi s pulsar with the K + producti on dFuchs et all 
120011: ISturm et al. L 1200 U lHartnack et all 120061) an d the tT/tt* 
production ratio measured in heavy-ion collisions (IXiao et all 
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2009), without resorting to exotic explanations such as a "fifth 
force' 'dWen et all l2009) or variations on the gravitation constant 
dWen et al.Ll2.012l) : it is enough to suppose that nucleonic matter 
undergoes a transition at high densities to a phase whose EoS is 
limited only by causality. 

On the other hand, comparing Tables|4]and[5]shows that the 
effect of quark deconfinement is to reduce the maximum possi- 
ble neutron-star mass, assuming that the speed of sound is lim- 
ited by cj V3. We stress, however, that the maximum density of 
neutron stars may not be high enough for perturbative QCD to 
be valid, and that the EoS of deconfined quarks might well be 
limited only by causality. 

It is interesting to compare our results with those corre- 
sponding to the assumption of a maximally stiff EoS, i.e., to 
the assumption tha t the EoS i s at the causal limit for all den- 
sities a bove »n dZeldovich l 1 19621: iNauenberg & Chapline . 
1971 iRhoades & RuffinT" 1 1 974t iMalone et al. L 1 1 975 ; 
Brecher & Caporaso 1 1 19761 iHegvi et al l 119751: jHartle L 1 1978 ; 
Lattimer et all Il990t iKalogera & Bavrri 1 1996c iKoranda et all 
1997; Sag ert et aU l2012). The results are given in parentheses in 
Tables [3] and |4] We found, however, that such configurations are 
thermodynamically unstable. Not surprisingly, the maximum 
mass thus obtained ~ 3.7M Q is very high. This upper bound has 
important consequences for identifying compact astrophysical 
sources. Taken at face value, it indicates that the soft X-ray 
transi ent GRO J0422+32, wh ose measured mass is 3.97 + 0.95 
M^ dGelino & Harrisonl 120031) . could be a neutron star. On the 
other hand, the requirement of thermodynamic stability imposes 
stringent constraints on the maximum neutron-star mass. In 
particular, the maximum mass is found to be reduced by about 
1.6 - 1.7M for the softest of our nucleonic EoSs, as compared 
to the maximally stiff EoS, as shown in Tableland Fig. [7] In 
this case, the identification of GRO J0422+32 as a neutron star 
is ruled out: it must be a black hole. 

If the core of neutron stars is made of deconfined quark mat- 
ter, the speed of sound in the maximally stiff EoS (obtained by 
setting «x = «n) will be limited by c/ V3. In this case, the in- 
equality .P q uai-k(«) > -PnucO) for n > « N , hence also £ qua i-k(«) > 
£nuc(«X is still satisfied but only in a restricted density range, 
since the speed of sound in nucleonic matter generally exceeds 
cl V3 at high enough densities. For this reason, the reduction 
of the neutron-star maximum mass after imposing thermody- 
namic stability is found to be much less dramatic, amounting 
to ~ O.6M at most, as shown in Tableland Fig. [8] 

As a matter of fact, the maximum masses of neutron stars 
without and w ith quark cores are well approximated by the se al 
ing relations (Hegvi et 



: cores are well approximated by tne sea 
ail 1 1 975t iHartle ll 1 978t IwTttenL IT984) 




9 10 11 12 13 14 15 16 17 18 
R [km] 

Fig. 7. Range of masses and radii of neutron stars with no quark 
cores (shaded areas), for maximally stiff EoSs and for EoSs 
satisfying thermodynamic stability. For comparison, we indi- 
cate the measured mass of PSR J1614- 2230 including esti- 
mated errors from Demorest et al. (2010). The dotted area de- 
limits the estimated r ange of masses of the soft X-ray transient 
GRO J0422+32 from lGelino & Harrisonl d2003l 




10 11 
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M max - 4.09 
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(21) 



Fig. 8. Range of masses and radii of neutron stars with quark 
cores (shaded areas), for maximally stiff quark EoSs and for 
EoSs satisfying thermodynamic stability. For comparison, we 
indicate the measur ed mass of PSR J1614 -2230 including es- 
timated errors from Demorest et all (12010b . 



(22) 



respectively, where Bq = 56 MeV fm 



5. Conclusions 



-3 



We have investigated the impact of a phase transition in dense 
matter on the structure of neutron stars, considering the stiffest 
possible EoS constrained by i) causality and ii) thermodynamic 
stability, i.e., the condition that at a given pressure the exotic 



phase should have a lower Gibbs free energy per baryon than the 
nucleonic phase. The latter condition is found to severely limit 
the maximum mass. 

Even if the phase transition is accompanied by a strong soft- 
ening of the EoS, we find that in the causal limit the maximum 
possible neutron-star mass is always increased above the value 
determined for a purely nucleonic EoS. In particular, the softest 
of our three EoSs, BSkl9, will then be able to support a neutron 
star as massive as PSR J1614-2230, provided the phase transi- 
tion begins at a density «n as low as 0.2 fm -3 . This shows, inci- 
dentally, that the existence of a two-solar mass neutron star is not 
necessarily incompatible with the soft nuclear-matter EoS that is 
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Table 3. Global structure of nonrotating neutron stars: maximum mass, corresponding radius, corresponding highest baryonic 
density of the nucleonic phase, corresponding lowest baryonic density of the exotic phase, and corresponding central baryonic 
density. 



Force 


M mm /M a 


R (km) 


n N (fnT J ) 


nx (for 3 ) 


"ccn (fnr J ) 


BSkl9 
BSk20 
BSk21 


2.03 (3.66) 
2.31 (3.64) 
2.42 (3.71) 


8.75 (15.62) 
9.99 (15.54) 
10.48 (15.80) 


0.20 
0.20 
0.20 


0.63 (0.20) 
0.49 (0.20) 
0.45 (0.20) 


1.41 (0.44) 
1.08 (0.45) 
0.99 (0.43) 



Notes. The quantities in parentheses refer to the corresponding results for the maximally stiff EoS obtained without imposing thermodynamic 
stability. 

Table 4. Same as Table[3]for hadronic-quark phase transition. 



Force 


M raax /M 


R(km) 


/?n (far 3 ) 


n x (fm- J ) 


np (fnT J ) 


n C en (fnr J ) 


BSkl9 


1.72 (2.26) 


9.95 (12.93) 


0.20 


0.34 (0.20) 


0.77 


1.24 (0.74) 


BSk20 


1.88 (2.26) 


10.90 (12.93) 


0.20 


0.29 (0.20) 


0.55 


1.17 (0.74) 


BSk21 


2.02 (2.31) 


11.70 (13.24) 


0.20 


0.25 (0.20) 


0.47 


0.90 (0.70) 



Table 5. Summary of results from Chamel e t al.l(l201 ll) for maximum neutron-star mass without an exotic phase. 



Force 


M ma JM B 


R(km) 




BSkl9 


1.86 


9.13 


1.45 


BSk20 


2.15 


10.6 


0.98 


BSk21 


2.28 


11.0 


0.98 



suggested by the measur ements of the kaon and pion productions 
in heavy-ion collisions dFuchs et alj | 200U ISturm et al. L [200 1 ; 
lHartnack et all [2006: Xia o et aUl2009l) . 

On the other hand, the presence of deconfined quarks in 
dense matter will generally lower the maximum mass of neu- 
tron stars, assuming the speed of sound is limited by c/ V3, 
as fou nd by perturbative QCD calculations at zero temper- 
ature (iKurkelaetall l20T(ih and non-pertubativ e lattic e QCD 
calculations a t finite temperatures (see e.g. iKarschl l2007t 
Bors anvi et all l20Tol and references therein). In this case, only 
the stiffest of our nucleonic EoSs, BSk21, will be consistent 
with the recently measured mass of PSR J1614-2230. If con- 
firmed, reported observations o f significantly more massive neu- 
tron stars with M > 2M n dClark et all 120021: iFreire et all 
2008; van Kerkwii k"etalll201 11) will hardly be compatible with 
the presence of qu a rk ma tter in neutron-star cores (see also 
lLattimer & Prakashl 1201 ll) unless the sound speed is signifi- 
cantly higher than cj V3. 

Considering the current knowledge of dense-nuclear matter 
properties, it would be difficult to understand the existence of 
neutron stars heavier than ~ 2.4 - 2.5M . This upper limit is 
considerably lower than estimates that did not impose the con- 
straint of thermodynamical stability. 
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